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Signatures of axion-like particles in the spectra of TeV gamma-ray sources 
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One interpretation of the unexplained signature observed in the PVLAS experiment invokes a 
new axion-like particle (ALP) with a two-photon vertex, allowing for photon-ALP oscillations in 
the presence of magnetic fields. In the range of masses and couplings suggested by PVLAS, the 
same effect would lead to a peculiar dimming of high-energy photon sources. For typical parameters 
of the turbulent magnetic field in the galaxy, the effect sets in at E 1 > 10 TeV, providing an ALP 
signature in the spectra of TeV gamma sources that can be probed with Cherenkov telescopes. A 
dedicated search will be strongly motivated if the ongoing photon regeneration experiments confirm 
the PVLAS particle interpretation. 

PACS numbers: 98.70.Rz, 14.80.Mz 
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I. INTRODUCTION 

One of the phenomenologically most important proper- 
ties of the hypothetical axions is their two-photon vertex 
that allows for axion-photon conversions in external elec- 
tric or magnetic fields 0,0. In particular, this coupling 
is used by the ADMX experiment to search for axion 
dark matter [3, Inland by the CAST experiment to search 
for solar axions 5|, . Generically such particles affect 
the propagation of photons in magnetic fields. For a lin- 
early polarized laser beam propagating in a transverse B- 
field, signatures are a rotation of the plane of polarization 
and the development of an elliptical polarization compo- 
nent 0, IS • The latter effect is also caused by the effec- 
tive four-photon interaction predicted by QED [l(| Q3[ • 

Recently the laser experiment PVLAS has reported 
such results with an amplitude about 10 4 times larger 
than expected from QED [12| . If one interprets this sig- 
nal in terms of photon-axion conversions, these measure- 
ments imply an axion mass m a w 1.3 meV and a cou- 
pling with photons g ai « 3 x 1CP 6 GeV -1 , where the 
coupling constant is defined in Eq. (JTJ) below. This com- 
bination of m a and g aj is incompatible with axions in the 
usual sense. Therefore, the new states require a differ- 
ent interpretation and are generically termed "axion-like 
particles" (ALPs), meaning bosons with a two-photon 
vertex where the mass and coupling strength are taken 
as independent parameters. 

The main problem with the PVLAS signature is that 
it violates simple astrophysical limits by a huge margin. 
ALPs are produced in the Sun and other stars by the 
Primakoff process where thermal photons convert in the 
fluctuating electric fields of the stellar plasma [l|, 0, • 
Assuming the PVLAS-inspired parameters, a standard 
solar model leads to an ALP luminosity so large that the 
Sun would burn out in 1000 years. Circumventing this 
vast discrepancy is the main theoretical challenge for the 
PVLAS particle interpretation [H, [H, EE El, OS HI • 

It is conceivable that the presence of the hot stellar 
plasma modifies the effective couplings or that these cou- 



plings are different at the momentum transfers relevant 
in stars. Therefore, it has been stressed that the PVLAS 
particle interpretation should be tested with experiments 
where the transition takes place in vacuum and where 
the momentum transfer is small 21]. Photon regenera- 
tion experiments ("shining light through a wall") are of 
particular interest because it will be fairly easy to con- 
firm PVLAS if the particle interpretation is indeed cor- 
rect. Several such efforts are now being discussed or are 
already under way [12, HE HI, notably ALPs at DESY, 
BMV at LULI, GammeV at Fermilab, LIPSS at Jefferson 
Laboratory, OSQAR at CERN, and PVLAS-regeneration 
at INFN Laboratory in Legnaro. 

If the PVLAS particle interpretation is confirmed, 
some radical new low-energy physics must be at work 
that prevents ALP emission from stars. However, other 
astrophysical settings provide conditions similar to the 
laboratory experiments, i.e., a vacuum environment and 
near-vanishing momentum transfers. One example is 
"shining light through the Sun" where a high-energy pho- 
ton source would become visible through the Sun by 
photon-ALP conversion in the solar magnetic field on 
the far side of the Sun, and their regeneration on our 
side . Another example is the double pulsar J0737- 
3039, where gamma rays emitted by one pulsar period- 
ically pass through the magnetosphere of the other on 
their way to us [261 ] . 

We here consider another example, the photon-ALP 
conversion in the turbulent magnetic field of our galaxy. 
Beyond energies of order 10 TeV, the gamma-ray flux 
would be depleted, leaving a distinct signature in the 
spectrum of TeV gamma-ray sources. Current data from 
Imaging Atmospheric Cherenkov Telescopes (IACTs) do 
not allow for a stringent constraint on this effect. How- 
ever, if the laboratory experiments confirm the exis- 
tence of ALPs with the properties suggested by PVLAS, 
this depletion must be included in the analysis of TeV 
gamma-ray sources by IACTs. Given the strong motiva- 
tion that would be provided by a positive laboratory ALP 
confirmation, dedicated efforts by present and future in- 
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struments would be mandatory that could provide an 
independent astrophysical signature of these novel parti- 
cles and/or allow one to study or constrain the turbulent 
galactic B field. 

We begin in Sec. [IT] with a summary of the formalism 
to describe photon- ALP conversions and turn in Sec. IIIII 
to phenomenological consequences on the propagation of 
TeV photons in our galaxy. In Sec. IIVI we briefly touch 
on the possible effect of millicharged particles on photon 
propagation. We conclude in Sec. El 



Equation ([3]) only holds in the perturbative regime where 
N Pq <C 1. For N sufficiently large, this result violates 
unitarity. It can be shown (Appendix of Ref. [28]) that 
the correct continuum limit after travelling over z 3> s is 



1 — exp 



2s 



(5) 



As physically expected, Eq. |5j) implies for z/s — > oo that 
the conversion probability saturates so that on average 
one third of all photons converts to axions. 



II. PHOTON-AXION CONVERSION 

Axion-like particles by definition have a two-photon 
coupling. For pseudoscalars, it is of the form 

L ai = -~g cn F iai F'" / a = ,g Q7 E • B a , (1) 



where a is the axion-like field with mass m a , the 
electromagnetic field-strength tensor, F^ v = he lxvpa F f " T 
its dual, and g ai the ALP-photon coupling with dimen- 
sion of inverse energy. For a scalar particle, the coupling 
is proportional to F^F^a. To be definite we limit our 
discussion to the pseudoscalar case, but similar conse- 
quences apply to scalars. 

As a consequence of this coupling, ALPs and pho- 
tons oscillate into each other in an external magnetic 
field. Under quite general assumptions, the probability 
for an unpolarized photon beam to convert to ALPs after 
traversing a magnetic field B = (B x , B y , B z ) from to z 
is (Appendix A) 



Pj — >a 
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dz'e- i2W / l °B x (x,y,z') 



dz'e- i2wz ' /lo B y (x,y,z') 



, (2) 



where for simplicity we have chosen the z-axis along the 
propagation direction. Further, Iq = ^irE /m 2 is the os- 
cillation length with m a the axion mass and E the photon 
energy. The meV range of ALP masses, relevant for the 
PVLAS particle interpretation, is so large that the pho- 
ton plasma mass is completely negligible by comparison, 
in contrast to the case of cosmic microwave conversion 
into intergalactic magnetic fields, studied in [27l. l28l] . 

We consider a simplified case where the field is of con- 
stant magnitude and random direction in each patchy 
domain, each with typical size s -C z, so that a large 
number N of domains is crossed. The previous expres- 
sion then further simplifies to (Appendix [A"]l 



P^ a {z) = NP , 
where the probability per single domain is 



g g 7 (|B| 2 )s 2 sinV S /Z ) 
4 (tts/Io) 2 



(3) 



(4) 



III. CONVERSIONS IN THE TURBULENT 
GALACTIC MAGNETIC FIELD 

For the PVLAS-inspired parameters m a — 1.3 meV 
and g a ~f = 3 x 10" 6 GeV~\ it is useful to write P in 
suitable numerical units. 



Pq = (1.556-B^GSpc) 
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Here, we have introduced g 6 = g ai /l0~ 6 GeV~ , B^ G 
is the root mean square (rms) magnetic field strength in 
micro- Gauss, s pc is the domain size in pc, m mc v is the 
ALP mass in meV, and E±q the photon energy in units 
of 10 TeV. In the second line we have replaced sin 2 with 
its average value \ because its argument is large and 
oscillates rapidly for any realistic energy resolution. 

Although the galactic B field has a regular compo- 
nent with several kpc coherence length, on small scales a 
turbulent component dominates (Ref. 29] and references 
therein). The power spectrum follows a Kolmogorov 
power law, with a lower cutoff at very small dissipative 
scales, perhaps as small as 6 x 10~ 4 pc [30| . but in any 
case at most comparable to 0.01 pc, with an rms inten- 
sity of order fj,G on pc scales. For typical galactic dis- 
tances of 10 kpc, there are approximately 10 6 domains 
with s « 0.01 pc towards a typical TeV gamma source 
such as the one at the galactic center [3l|, HH, HH, [34[ . For 
nominal values of the parameters in Eq. ([6]), Pq ~ 10~ 7 - 
10~ 6 at energies of order 10 TeV, implying N Pq ~ 0.1-1. 
Therefore, observable effects must be expected. 

In Fig. Q] we show the spectral modification of a TeV 
source at the galactic center (distance 8.5 kpc) superim- 
posed with H.E.S.S. data [35|]. For illustration we have 
used Eq. © with g 6 = 3, B^ G = 0.7, s pc = 0.01 and 
m mc v = 1 and we have assumed that the power-law 
spectrum does not break before 60 TeV. Photon-ALP 
oscillations (dashed curve) cause a downward shift of the 
spectrum at high energies, i.e., a change of normalization 
of the typical power-law spectrum (continuous curve) be- 
tween low and high energies (E > 10 TeV). The maxi- 
mum shift is 33% when the conversion saturates. 

Evidently current data do not allow for a serious con- 
straint on this depletion effect. Note, however, that 
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FIG. 1: Spectral energy density E 2 x dN/dE of photons from 
the galactic center source, for the 2004 data (full points) and 
2003 data (open points) of H.E.S.S. [3f|. Error bars represent 
95% CL. The continuous line shows the best-fit power-law 
dN/dE ~ E~ T with Y = 2.25 @. The dashed line shows 
the effect of photon-ALP conversion with coupling and mass 
suggested by PVLAS. 



the large error bars at high energy are only due to 
a lack of statistics. The points reported in Fig. Q] 
are based on 17 hours of data in 2003 with two tele- 
scopes and 48.7 hours in 2004 in the four-telescope ar- 
ray mode. Already the current generation of IACTs 
(CANGAROO-1TI, H.E.S.S., MAGIC, VERITAS) may 
have a sufficient aperture to probe this scenario, if dedi- 
cated campaigns were motivated by a positive laboratory 
detection. 

To be more quantitative, one would model the turbu- 
lent field as a Gaussian random field with zero mean and 
an rms value B Tms [2!| . Each of its components can thus 
be written in terms of its Fourier transform as 



(27T) 



(7) 



where the phases <fo(k) are random. For an isotropic and 
homogeneous turbulence, the Fourier modes satisfy 



B 2 (k) 
8irk 2 



^)-2;,-!'^(k k'i 



,(8) 

where the tensor in brackets implements the condition 
V • B = 0. In the generic case of a power-law spec- 
trum with index a between the scales s m ; n and s max , i.e., 
between wavenumbers = 27r/s max and fcn = 27r/s m i n , 
one has 



B 2 (k)=B 2 ms (a~l)k- a (k 



1-a 
L 



(9) 



which is already normalized such that (|B(x)| 2 ) 
In the limit <C fcH, and if a > 1, one finds 

B 2 (k)^B 2 ms (a-l)fc-X" 1 - 



= B'i 



(10) 



Therefore, the field averaged over scales less than s is 

(|B(x)| 2 ) s = ^ ms(s / Smax )-i. (n) 

For the Kolmogorov spectrum a = 5/3 suggested by the 
data, this means that the rms intensity of the field varies 
as s 1 / 3 . The intensity below 0.001 pc is then only a fac- 
tor 10 weaker than the //G level at the pc scale. Below 
0.01 pc it is only a factor ~ 4 lower than at the pc scale. 
Therefore, our simple estimate of Pq may be too opti- 
mistic by an order of magnitude. However, the effect 
would still be observable simply by looking at a factor 
~ 3 larger energies. Additionally, the true field config- 
uration may be more complicated, and recently a more 
intense turbulence than previously estimated has been 
suggested [36j . 

In a more detailed treatment one would consider 
stochastic realizations of the realistic power spectrum of 
the turbulent B field. However, for our purpose sim- 
ple estimates are probably more instructive and show 
that: (i) Possible effects may start manifesting them- 
selves around 10 TeV, and are more and more likely to 
show up at 20-30 TeV. (ii) The smaller the characteris- 
tic scale of turbulence of a given intensity, the larger the 
number of domains available, and the lower the energy 
at which the effect appears, (iii) The conversion prob- 
ability depends on E 2 . Therefore, on the scale of the 
typical energy resolution of a Cherenkov telescope, the 
depletion rapidly drops from negligible to the saturation 
value of 1/3. (iv) The phenomenology described here 
would be universal, affecting both galactic and extra- 
galactic sources. Yet, the exact energy at which the shift 
manifests depends on the properties of the field along 
that line of sight. Although for all sources the light must 
cross the galactic B field to reach us, one may not ex- 
clude an additional role of a small-scale field close to the 
sources. Our estimate for the onset of the effect is conser- 
vative, especially for extragalactic sources. As a general 
rule, for sources in similar directions, the more distant 
ones may manifest the signature at lower energies. 



IV. MILLICHARGED PARTICLES 

Another particle-physics explanation of the PVLAS 
anomaly postulates the existence of low-mass milli- 
charged particles [13, HI]. We briefly check if this hy- 
pothesis would also affect the propagation of photons in 
the astrophysical context. 

At TeV energies, the extragalactic medium becomes 
opaque due to the onset of e pair-production on the 
diffuse low-energy photon backgrounds. At a few PeV, 
the mean free path of photons reaches a minimum of 
X e % 10 kpc due to pair production on the Cosmic Mi- 
crowave Background (CMB) [39]. The threshold energy 
Eth ~ 3 x 10 14 eV scales as m 2 and the cross section as 
e /wig. Scaling these quantities to millicharged particles 
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with charge g«e and mass m q <C m e one finds 




and 

^= E -(S) 2 - (13) 

The preferred mass range of the millicharged candidate 
is 0.01-0.1 eV, i.e., m q /m e - 2 x 10" 8 -2 x 10" 7 . The 
peak of the cross section is very close to the threshold and 
would fall in the (10 -16 -10~ 14 ) x 10 15 eV range, i.e., rang- 
ing from infrared to ultraviolet. Sources at cosmological 
distances do not show such a universal dimming. The 
conservative requirement X q > 1 Gpc implies q < 10~ 5 e. 

A much more constraining limit of q < 10~ 7 e arises 
from spectral distortion effects of the CMB that may al- 
ready rule out the millicharged particle explanation of 
PVLAS [13]. In any event, it appears safe to assume 
that millicharged particles with the relevant properties 
would not affect TeV photon observations. 



V. CONCLUSIONS 

The unexpected optical properties of the vacuum sug- 
gested by the PVLAS experiment has inspired various in- 
terpretations in terms of axion-like particles. The severe 
conflict with stellar structure arguments implies that this 
interpretation requires radical new physics at low ener- 
gies. If the new particles interact differently in a stellar 
plasma or at vanishing momentum transfers, they may 
still show up in the upcoming photon regeneration ex- 
periments. In this case one necessarily expects signatures 
also in other settings that are characterized by a vacuum 
environment and/or small momentum transfers. 

We have discussed possible signatures of PVLAS par- 
ticles in the spectra of TeV gamma-ray sources in our 
galaxy. If the PVLAS signal can be attributed to photon- 
ALP conversion in the laboratory, the same effect must 
occur in the astrophysical context. For an ALP mass 
around 1 meV, as suggested by PVLAS, one would 
observe a peculiar distortion in the photon spectra at 
Ery > 10 TeV due to conversions in the turbulent galac- 
tic B-ficld. This process would take place under better 
vacuum conditions than are achievable in the laboratory 
and the momentum transfer would be extremely small. 

Present data from TeV gamma-ray telescopes do not 
allow for a stringent constraint on this effect. However, 
a positive ALP detection would strongly motivate a ded- 
icated search, perhaps allowing one to find signatures for 
ALPs in current or future instruments and to investigate 
or constrain the properties of the turbulent magnetic field 
in the galaxy and beyond. 
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APPENDIX A: PHOTON-ALP CONVERSION IN 
A RANDOM MAGNETIC FIELD 

We here derive the photon-ALP conversion probabil- 
ities [Eqs. (HI)-©] in a random magnetic field distribu- 
tion. These detailed results are not used for the simple 
estimates derived in our paper, but would be necessary 
for a detailed treatment involving the numerical study of 
different realizations of the turbulent galactic B-ficld. 

For relativistic ALPs, the equations of motion follow- 
ing from Eq. (fTJ) reduce to the linearized system Q 

{oj-\d z + M) {^A y j =0, (Al) 

where z is the direction of propagation, A x and A y are or- 
thogonal components of the photon field in a fixed frame 
perpendicular to z, and to is the photon energy. The 
mixing matrix is 

(^xx ^xy ydwy-^x \ 

i \x i A yy \g ai B y , (A2) 
^day^x "29a~{By A a J 

where A a = —m 2 a /2Lo. Notice that the component of 
B parallel to the direction of motion does not induce 
photon-axion mixing, since only B x and B y enter the 
third row/column of M.. The entries = x, y) that 

mix the photon polarization states are energy-dependent 
terms determined by the properties of the medium and 
the QED vacuum polarization effect. We will neglect the 
latter because it is sub-dominant here. 

The Ay terms have a simple interpretation when the x 
or y direction coincides with the transverse field direction 
B T = B — (B • z) z. We can then specify the previous 
equations for the case of a single domain with uniform 
magnetic field Bt, whose modulus will be denoted by 
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Bt = |Bt|. Equation (IA1|) is in the new basis 



A± 

(uj-id z + M) ( Ay 
a 



= 0, 



(A3) 



where i = _L or 
matrix is now 



refer to the Bt direction. The mixing 



(A4) 
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All 


A a7 





A Q7 


A a 



where 



A„ 



A p i 
A p i 



CM 



A_l" 



A p i = 



2 9 ay 



B, 



(A5) 



Here, a-i 2 j = Ai:an e /m e is the plasma frequency with m e 
the electron mass and a the fine-structure constant. The 
Faraday rotation term Ar, which depends on the energy 
and the longitudinal component B z , would couple the 
modes Am and A±. While it is important when analyzing 
polarized photon sources, it plays a negligible role here. 
The Acm terms describe the Cotton-Mouton effect, i.e., 
the birefringence of fluids in presence of a longitudinal 
magnetic field, with \A9 M - A5 M | oc B\. These terms 
are of little importance for the following arguments and 
will be neglected hereafter. 

Therefore, we finally concentrate on the simple two- 
level mixing problem 



- i<9 z 



A p i A an 
A a7 A a 



An 



II =0. 



(A6) 



The solution of this system follows from a diagonalization 
of the mixing matrix by a rotation with an angle 



i) 



arctan 



2A ai 



A p i - A a 



(A7) 



In analogy to the neutrino case, the probability for a 
photon emitted in the state A\\ to convert to an ALP after 
traveling a distance s in a constant transverse magnetic 
field Bt is 



Poh^a) = \(A {[ (0)\a(s))\' 

= sin 2 (2tf)sin 2 (A osc s/2) 
^2 sin 2 (A osc s/2) 
" (Aa ^ Sj (A oscS /2) 2 

where the oscillation wave number is given by 

AL = (Api-A a ) 2 + 4A 2 . 



(A8) 
(A9) 

(A10) 



(All) 



The conversion probability is energy independent when 
2|A a7 | » | A p i — A a | or in any case when the oscilla- 
tory term sin 2 cc/a; 2 « 1 in Eq. (| A10|l . corresponding to 
A osc s/2 < 1. 



We now return to the 3x3 formalism to derive a per- 
turbative solution. In a fixed x-y-z frame with z the 
direction of motion, the propagation equations are 



= 0, 



(A12) 

where c 7 = cos 7 and s 7 = sin 7 with 7 the angle between 
Bt and the y axes (measured clockwise). Further, from 
Eq. (A5) one can write 
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(A13) 



The field strength entering A a7 is Bt = |Bt| = 
|B| I sin^|, where ip is the angle between the field and the 
photon propagation direction. Thus we have B x — Btc 7 , 
By = BtSj that are all z-dependent quantities. All of 
the Aij are z-dependent as well because this applies to 



7, n e , and Bt, entering the quantities in Eq. (|A5p . 

Since the ALP is weakly coupled, the 3rd row/column 
off-diagonal terms are much smaller than uj, and it makes 
sense to write 



i<9 2 A = (Ho + Hi) A 
where A = (Ax, Ay, a), 



and 



Hi 



H =u)I 



A p i 
A p i 
A a 



(A14) 



(A15) 



^0,7 ^7 ^07 ^7 



(A16) 



For gay — * this equation is solved exactly by A^ ) (z) 
Uo(z)A(0), where 



Uq (z) = exp 



-i / dz'Ho(z') 
Jo 



(A17) 



If we now include the perturbation, the complete so- 
lution can be written perturbatively in the interaction 
representation. In particular, to first order we have 



A int = U ] A, H^t = UIHiUq, and 

A$(z) = -i f dz'H- mt (z')A^), 
Jo 

and Afa[.(z) = A(0) because 

A^(z)=U^(z)=U^U A(0). 



(A18) 



(A19) 



Since Hq is diagonal, Uq has the general form Uq(z) 
diag[e~ iQ ( z ),e~ ifc ( z \e- ic W] so that 
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The ALP amplitude developed at distance z is then 



e" i ( c - a )A a7 s 7 

H^t = ( e-^-^A^Oy 



^ ' ^ ^^Q'ySry ©' ^ (-J >-y C<-y 



(A20) 



a (D( z ) = _i^l / «I; 



'{A x (0)B :i; (z , )e i[c(2 ' ) - a ^' )1 + A y (0)B a (z')e i[c(z ' ) -'' (z ' )1 } . 



(A21) 



This result is a straightforward generalization of the one derived in Ref . [9|] . The probability for photon- ALP conversion 
is then schematically 



P^ a (z) = lA^HZil 2 + \A y (0)\ A \1 2 \ 2 + 2Re[A x (0)A y (0)l 1 T 



(A22) 



For an unpolarized source, an average over the initial state has to be performed. The interference term averages to 
zero, and (|^(0)| 2 ) = (|A,(0)| 2 ) = 1/2. Then 



or equivalently 



(A23) 



P-y^a(z) 



gg 7 |B| 2 



dz'sinV'(2 , )e i(A »- A '" )2 'c 7 (z') 



dz'sinV'(z , )e i(A °- A P l)z 's 7 (z') 



(A24) 



r 



Here we have assumed A p i to be independent of z. 

We next consider a "patchy" pattern of domains of 
equal size s and constant field in each of them. We will 
show that, when evaluated after a distance z w N s, with 
N >■ 1, the conversion probability is roughly the product 
of the conversion probability in a single domain times the 
number of domains. Except for the replacement s 7 — > 
c 7 , each one of the two integrals in Eq. (|A24p can be 
evaluated as follows, where lo = 2tt/(A p \ — A a ), 



1 = 



dz'sinV>(^)e~ 27ri2 ' /io s 7 (^) 



k=i •' Zk 



• 2 f nS 
* 8111 l "fc 



(2z fe +s)/*o 



fe=l 



interference 
terms 



(A25) 



Here, N = z/s ^> 1 and /ifc = | sin^fel s 7 (&) or /ifc = 
| sinV'fcl c 7 (fc) is a random variable in the interval [—1,1]. 
The random nature of the field directions implies that 
the interference term vanishes on average. For the ge- 



ometrical factor we have (/i 2 ) = (sin 2 tp sin 2 7) = 1/3. 
Then we find 



P^a(z) 



5 2 7 |B| 2 ll . 2 fns 



8 7T 2 



sin 



„ N 
, x 2 x — 

lo J 3 



N ((A Q7 )s) 



2 sin 2 (|A pl - A„|s/2) 



(|A p i-A |a/2)2 



(A26) 



having the structure of a probability per single domain 
Po times the number of domains N. We stress that 
Eq. (|A26[) only holds perturbatively, i.e., (A aj )s <C 1 
is a necessary condition. 

In the limit |A p i - A tt | > (A a7 ), we have in Eq. (|A10j) 
that A osc = |A p i — A a | and Eq. (|A10|) coincides with 
Eq. (jM6|) . provided that B T -> (|B|) = |B|/V3 be- 
cause of the projection effect. In the opposite limit 
|A p i - A„| < (A a7 ), Eq. IpggjP reduces to 



P ( ^ a (z) » N((A a7 )sf 



(A27) 



again in agreement with the corresponding limit of 
Eq. (TATOl) . 

This exercise shows explicitly how the classical rule of 
"adding the probabilities" instead of amplitudes arises 
from the randomness of the polarization and of the field 
configuration over scales much larger than s. However, 
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since we used first-order perturbation theory, the valid- 
ity of these results breaks down when P 1 ^ a {z) becomes 
large. This is always the case for z large enough, since 
we are not including the back reaction a — ► 7, that are 



second order in g ai and that prevent the violation of 
unitarity. In the saturation regime, the correct general- 
ization of Eq. (|A10|) is provided by Eq. ([5]) as discussed 
in the text. 
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